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1. INTRODUCTION 

At present large number of electronics devices based on heterostructures. Heterostructures could be 
grown by using different well-known approaches [1-12]. At the same time at large number of experimental 
studies on growth of heterostructures [1-12] we find small quantity of theoretical works with prognosis of 
processes [12]. Main aim of this paper is to analysis changing of properties of epitaxial layers with a 
changing of the values of different parameters of the epitaxy process with account chemical interaction 
natural convection. 

We consider a vertical reactor for gas phase epitaxy (see Figure 1). The reactor includes into itself 
(i) external casing; (ii) keeper of substrate with a substrate; (iii) spiral around the shell in the region. The 
spiral generats an induction heating to activate chemical reactions. The heating activate chemical reactions 
between materials during epitaxy. A mixture of reagents in gas phase enter to inlet of the reaction chamber in 
atmosphere of a gas-carrier. Our aim framework this paper is analysis of changing growth process with 
changing of physical parameters and parameters of technological process. An accompany aim of our paper is 
accounting of natural convection and chemical reaction. The third aim of our paper is development of 
analytical approach for prognosis of growth process to increase adequateness of the prognosis. It should be 
noted, that analytical approach for modeling are usually more demonstrative in comparison with numerical 
one and leads to decrease calculation time. 
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Figure 1 (a). Structure of the considered reactor for epitaxy Figure 1 (b). View of keeper of substrate 
from gas phase with sloping keeper of substrate from side with approximation of the keeper 
by sloping lines with angle of sloping qi 


2. METHOD OF SOLUTION 

We solve our aims by determination and analysis of spatio-temporal distribution of temperature and 
concentrations of reagents in the considered reactor. The spatio-temporal variation of temperature was 
obtained by induction heating due to electric current in spiral around keeper of substrate with a substrate to 
activate of chemical reactions gas-reagents before substrate during formation of required epitaxial layer on 
the substrate. We determined the considered temperature by solution of the following boundary value 
problem [13]. 


orn ba) _ div(A : grad[T (r, $, z, t)] — [?QCr, $,zt)- u(r, $,z, t)] -c(T) -T(r, bz, t): 
Cs(r, $, Z, t)} tp, $,z, t), (1) 


Here is the flow velocity of the gases-reagents; c is the heat capacity of the system; T (r,,z,f) is the 
distribution of temperature in space and time; p (r,z,t) is the heat power density in the considered reactor; r, 
9, z and t are the current cylindrical coordinates and time; C»(r,,z,f) is the distribution of the concentration of 
the mixture of reagents in space and time (we assume, that two reagents with concentrations of inlet in 
reactor with concentrations Ci(r,,z,f) and C» (r,,z,t)); C3(7,,z,t) is the material of the epitaxial layer in the gas 
phase; A is the thermal conductivity coefficient, value of this coefficient has been determined by following 


relation: 2=Vic,p/3, here v - J2kT/m describes modulus of the mean-square gas molecules velocity; / 


describes mean free path of gas molecules between collisions, c, is the heat capacity of the gas at constant 
volume, pis the gas density. 

To solve above equations we shall to take into account moving and quantity of reagents. In this 
situation we shall to solve equation of Navier-Stokes with the second equation of Fick. We assume that the 
radius of the substrate holder R larger, than diffusion and boundary layers thickness. We also assume that the 
flow of gas is laminar. In this case the above equations could be written in the form 


ar + GV) = -v (2) + vas, (2) 

ane = div{D, : grad[C,(r, ¢,z,t)] —[00,¢.2z,0) - 9] Gr, 62,0] - là oz O- 
6 (r, $, Z, t), (3a) 
ee = div{D, - grad[C,(r, ¢,z, t)] - [7 ¢,z,t) - 9| -C,(7,¢,2z,0)} — km 2,0) - 
C; (r, $, Z, t), (3d) 
Tun = div(D, - grad[C3(r, ¢,z, t)] - [(7 p, z, t) — 9| C8 6,0] ks oz t) + 
C3 (r, $, Z, t), (3c) 
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Here D; are the coefficients of diffusion of reagents of the considered reagents. As an example of 
reagents, we consider trimethylgallium (CH3)3Ga and arsenic hydride AsH3. As result of reaction we consider 
GaAs. As gas-carrier we consider hydrogen. P describes pressure of gas in the reactor; v describes kinematic 
viscosity of reagents. We assume, that all molecules will deposit on substance. Based on this assumption we 
used the following boundary and initial conditions. We also assume, that reactor is cylindrical and d»««di. In 
this situation we can write 


Ci (r,9,-L,t) = C2 (7,9,-L,) = Co, C3 (r,9,-L,1) = 0; Ci (r,- 9uz,t) =C (r Quz.) = Ci (r,z- Q1,Z,t) -C 
(rz *9uz); Ci (nez0) = C2 (7,9,2,0) = Cod (Z+L), Cs (ro,20) =0; Ci (0,920) + o; 
Ci(r dz, t)ls = 0; T (r.- uz) = T (9uzt) = T (rz guzt) = T (r, zot); T (92,0) = T, 














OT (r,,z,t) 4 OT (r,¢,z,t) OT (r,,z,t) Ovr(1,¢,z,t) 
-1 ==] -oT*(R,$,z t); = ; — | = 0; 
| Rp ) ag $-0 ag p=2n or r=R 
POOPED! = 0; ATOR! = oTr, p, -L t); POP! = 0; v (pun = vr (Pz) 
or T-R Oz S or r=0 








= v, (rz- Quz) = vr (rzt gut); T (0,9,2,0) # ©; vo (r.-Quz.t) = Vo (7,01, z.f) = Vo (1,4-G1,2,t) = Vo 
(r,2+1,2,t); Vz (rs quz,t) = Vz (7,Q1,2,t) = vz (rz-quzt) = vz (rum euzt); vr (1,9,-L,0 = 0; vr, 9,L,0 = 
0; v. (0,@,2,0) # oo; v(r,Q,-L,0) = Vo; vrtd2/2,9,z€[-d2/2,d2/2],0) = @z-cos (W)-tg (915 vo (r.p, L.t) = 
0; vo(0,Q,z,t) # ©; v{r,¢—,0,1) = 0; vr, o, Lt) =Vo, v{r,9,L,t) = Vo, v{0,9,z,t) # oo; v(r,9,z,0) = 0 
ve(r,Q,z,0) =0; vdr, 9,0, = or, (4) 


Here o =5,6710° Wm?K*, T, describes room temperature, œ describes frequency of rotation of 
substrate. Using cylindrical system of coordinate leads to the following form of equations for projections of 
reagents velocity 

















dv, 18 [ dv-(7,¢,z,t) ES eg ae) a2 v,(r,o,z,t) av, Vg AVE Ov; O0 fP 

We y (E trant), A orgao reae] y de YAM y 2 (P) (sy 
ôt r or 09? 0z? or r 0$ Oz or Vp 

ðvr _ 1 à | meern y 1 E Nas a? ve lr, pz t) dvr VeVe ðvz 1 ð (5 
at vf ar |” ap? az? "ror r 8$ — az p (OP) 
Ov; _ (1 A [.OvzG, 5920] y - 1 E STD) Tutgan) Avr VeVe ðvz ð (5 

at — v ol 8$? az2 ar rs 8$ Vida az pi (Se) 


Now we will calculate solution of (5) by method of averaging functional corrections [14-19]. 
Equations for first-order approximations v,1, vio, Vz1 of the considered components takes the form 


vir _ a E) avg __ 12 C) viz _ a E) (6) 
a | OrNoP 8t — rəææ\p? at OzNo/ 


Now we integrate both sides of the above equations. The result of integration could be written as 





ð tP 
-dT, Vig = 37 0p 


Vir = -dt, Vig = 


tP 
^ ar zs p Saas ale p (7) 


Equations of the second-order approximations of the velocity projections vm, vig, vz1 could be 
written in the form by using standard procedure of method of averaging functional corrections 











Over _ [ a ( arar) 1 0747 aa a E) ova (a2g+V1g) Avir 
ae par d ar Tz ag? i az2 ar \p (azr + MG r ap 


Par, 





(a2; t Mids 


(8a) 








avrg __ [103 ( rie) 1014 , 07r14 18 (5 and (225 *v19) Vip 
at -—PEE Pe) tz ag? az? ap Np (azr + vir) r ag 


0 
(azz viz) =, (8b) 

















üv2z [2 a ( araz) 1 024, add a E- L (a2g+01g) viz 
=y r = —— a, T vir) — - ———— 
at rar ar r2 a? i az? az Np (Gar it) g 

mt 





(a2; + Viz) 


(8c) 
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Integration of the above equations on time ¢ leads to the following result 

















Ovir 1 avar 0?vi,y t Ovir 
Yh IG x)te age wae" "ES dr) — fy (aas + Mir) ge 
t pes Ovir t Vir 
d i 8$ dt — d: (azz + Viz) az dt, (8d) 














(r=) 10?v,g Prip 18 tP t 
Yap — V f, E d )fz 32 aga dr — (opde) - fo (tae 


t (@oagtv Ov Ov. 
fy aeria) ab de filar + vi) WP dr, (8e) 


r 











Qviz 10?v,2 , O0?vgz a tP t Oviz 
2; 7 Vo tv f, Fst bn ag? T ae] d 5; (o t) Jo @ar + Vir) Ge at — 


v +019) viz avy 
Spe edet] ae vu edu. (8) 


Average values œ, 029, 0»; have been calculated by the following standard relations 


2 O KR 2 L 
a» = — Hh [s £p E Var — Yır)dzdġdrdt, = $7 sand if r f, "Sze E 
v4 )dzdódrdt, a, = — f" f rfe f Ozz- aes (9) 


cT 


Here O is the continuance of technological process. Using of the first- and the second-order approximations 
of the velocity projections in the above relation leads to necessity of solution of the algebraic system of 
equations, which is presented bellow 


41025, + By dog + C105; = D, 
4202, + B305$ + C22, = D; (10) 
A305, + B305$ + C505; = D3 


where, 








A 21*fe(0-Of rf, f, dzdodrdt, B, = f, (0 —O fy f; f. e  dzdódrdt, C, = 









































L ab 
2m OvVir 1 0^vir j genie 
C, 258 RD, =v x [PO -O fr fof, Pa (r ER MM | azdedrat — 
fj (0 Of rf; f. vi, SE dadpdrdt — = ?R?v$ — f; (0 — €) f, TI Vio ot dzdddrat, 
A; f, (0 - O fo rfe. f, zt dzdpdrdt, B, 2 1 f, (0 — t) f, id de; oe da dpdrdt, 
2 ðv 1 0?v a*v 
D,-vff(6- Off rf?" f Eç rey ME re | debra gen 
2 avar 2 avar 
DS rN Si Mir St dzdédrdt -Z0 RV — (6-0 fo fp f vie Ta dzdpdrdt, 
cae TORV, As = fo (O-—t) fir fo” farm dzdódrdt, 
- f? (8 -Of rx f. P dz dódrdt, 

E 2T Oviz 1 Swan y a? s 2y2 _ 

ur OKT he Ex 2M x pedes 9? R? Vj 








fj (6-0 Sr" SÉ vi s dzd$drat — f (0 — t) ff [E [uds p dzddrdt. 
Solution of the above equations could be obtained by standard approaches [20] and written in the form 

azy = A,/A, Ong = Ag/A, 055 = 4/4, (11) 
where 

A = A,(B,C3 — B3C,) — By (A563 — A3C2) + C, (A5Bs — A;B;), 

Ar = D,(B;C3 — B3 C2) — Bı (D2C3 — D4C2) + C,(D B4 —D; B;), 


Ag = D,(B2C3 — B3C2) — Bı (D2C3 — D3C2) + C,(D,B3 — D3B2), 
A, = A;(B,D3 — B3D,) — Bı (A2D3 — A3D,) + D,(A2B3 — A3Bz). 
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In this section, we calculate the required projections of the considered velocity by the method of 
method of averaging functional corrections. Usually, the such approximation is enough to obtain main results 


of analysis [15-19]. 


Now we will analyze distributions of temperature and concentrations of components of gas-reagents 
in space and time. In this section we will consider (1) and (3) by using cylindrical system coordinate 


aT(r.z,t)  4,08?T(r$zt) 1 0?T(r,¢,z,t) PT = 
Coa À z tas 39? +2 + p(r,¢,z,t)—c- 





2 (ur, $,z, t) — 


v.(mr,ód,zt)l-Cmézto0-Tmozt- £2 (DC. $,z, t) — vs(r.d,z, 5] -Cs(r,¢,Z,t) - 


a - 

T(r,$,.z t)) - c- s; vs gz t) — 9,06, zt] Clr, pz O TG pz 0). (12) 

óCci(rdzt) 10 aC(rpz.t)] , 1 9 [py A(z.) " 

at ^ rór |, or T r2 op A 0o ^ rÓr sra. $, 7 e [v.(r, $, Žr D 

ac t 

v?.(r,$.zt)h--— = [p Es 23 2 (rc, zt): [ver $,z, — 9. pz t} - 
kr, ó, zt), é zt)- 2G, $,z,t)-[v,(r, , z t) — (7, bz, t)]}, (13a) 

üCcQmzt) 10 0C (r,biz.t) 0C r.biz,t) 

ein Li [po teana zz = |p, aaa -+2 (rey. d, t): [v.(r, $, z, t) — 


Br pz ON + 2 [Dp 85929] - 17 (reo Gr, 0 - [var 6 2,0 - dor. $,2,0)]} - 


à : 

k; (r, $, Z, t) t Cyr, $, Z, t) p az (C; (r, $, Z, t) : [v Gr, $, Z, t) "S v, $, Z, t)] (135) 

acro zt) 18 acs(r.pzt)] , 1 3 [p 3Cr ezt) _ 
ðt — ^ rór EP ar RS 9e Ds ap ^ rór Zirc (6,20): [gz t) 


(r, ġ,z,t)] JE |p; ar | ta res. pz t): [vg d, z 0 — 94,0, 2, 0]] - 


k3 (r, $, Z, t) " C3(r, $, Z, t) M < {C3 (r, $, Z, t) 9 [v (Cr, $, Z, t) ux vr, $, Z, t)]}. 


(13c) 


To calculate distributions of temperature and concentration of reagents in space and time, we use the 
same method of averaging functional corrections as for solution of (2). We used the same algorithm of 
method of averaging functional corrections to solve (12) and (13). The result of calculation of the first- order 


approximations 7'(r,,z,f) and Cii(r,,z,t) of considered distributions could be written as 


T, (r$, z, t) =T, f; C972 qr — arac jesus cuneo anl pcs 

















or 
a s 
arac ft [ve (r.o. eet iD, P dt = gu fi al ele ban terr dt, (14) 
Ca, o, zt) = cuc i Poir rb ten dt — ay f ORO Le C, — 
ð 9,21) -ve (rez, 
euc (t Peto rainten] desea ge LE aey, dı (15a) 
Cor, d, z, t) = — 92€ fe Topo s + Co — G21c (1 $,z, v)dt — 
ð Q, 
a [i persana] g ae n? UP -PTET dr, (15b) 
Ca (r,ġ,z t) = Tac f neo eeN qu + Q31¢ jn k3(r,$, z, T)dt — 
a durdjepl dum Lie 
Eno Pe pee een dt— ts, ie Plea ren) td do (150) 


where dir and ajc are the average values of the considered approximations. Now let us calculate these 


average values by using the following relations 


1 O pR 2 L 

a =— zh hr f, STC, ¢.2z,1)dzdpdrat, 
1 O pR 2 L 

Gic=—a lS rh " [Cal Q, 2, Ddzdpdrdt 


(16) 
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We calculate average values air and aiic by using the following standard relations 


Bic = Co/L- {1 + —— (0 — 0f A [v,(R, $, z, t) — 0, (R, d, z, t)]dzdódt + m, 
Qa = [r. t—À0-0frf p ( + S (0 — 
Of flv. o, zT) - 90, z ]dzdódt -—-— 0-A N fp fv oz t) - 


1 
8, (r, d, z, t)ldzdgdrdt ] — x (0 tdt — 2% 4a) ) 








TORL 


We calculate approximations of temperature and concentrations of gases with the second-order by 
method of averaging of function corrections [14-19]. The result of calculation could be written as 


cT, dz t) =A fL 0979 qp c. P (o, p, zT) 6,7 p, zT): laze + Clr, o.z, 0] - 





[a27 + T, (mr, $, z, )])dv T e dr h ([aac + Ca (r, $,2,7)] - [vor Q, 2,7) — 
elr, p,z,T)| - [aor + Ti (7, Q, z, 1)]}dt +T, + (16) 





44 gt O92 qr Ce 9 = f vs, $,z,1) — (r, 6, 2,T)] + [dae + C (7, o, z, 7)] + [mor + 
a $,z, T)]}dTt tf As $ z, t)dt, C (rd, zt) = Lf rp 95929 dr 4 
7 aah D aan dt + az fo D PC Cp) dt + Coó(z + L) — —a12c f ki(r, 6,2, v) [aac + 


Car, $, Z, ids m Lr f [aise + CQ, $.z, 1)] t [v,(r, $, Z, T) E» v.(r, $, Z, t)|dr} zr (17a) 





Soo laze + Ci o, zT) [vo Gr. $,z,t) = vor, $,z, 7)]d: -2 filarze + Clr, $,2,7)] - 


[P(r 2,2) — 9, (m. sz, Ddr, Clr, ds zt) = 1. fy DOPED de T prp eodem dr 4 








$ D Pa C022 dr + Co (z + L) — a22c f kar, $,z, v) [a22c + Cai (7, 4,2, 1)]d1 — 
73 t e 
= sr rÍ [@22¢ + Cat, $,z, 1)] Ë [v.(r, $, Z, T) E v.(r, $,z,c)ldr) 7 (17b) 


rór 


e [G22c + C316, 2,2] < [vo (r, p, z, T) — 94m, 2,7) |dt — az fle + Coi(7,,2Z,T)] - 


[v,(r, $, z,T) — v,(r, p, Z,T)|dt,C32(r, $, z t) = is A rp 2630922 qt 303 Fog! D -—— t 
à rt Cai (ribuz) 
zh poU dr + Q32c fe k3(r, Q, z, PM + Car, $, z, 7) dt — M rf [@32c + 


Oz 
C34, $, Z, 1)] t [v (r, $, Z, T) x v.(r, $,z, 7)Jdr] Te Sa fy sac + Car, $,z,1)] $ [ve Gr, $, Z, T) E 
Var, $,z, t)|dt — 2. fo [ta2c + C31(7,¢,2,7)] - [v (m $, z 7) — v;(r, , z, 7) ] dt. (17c) 


Average values of the above approximations were determined by standard relations 


1 O pR 2 L 
ar — us hr h " [02 — T0dzdéárát, 


1 0 cR 2 L 
Ore =a fy hr h S- (C2 Codzdédrat. (18) 


We calculate average values oor and aizc by using the following standard relations 





"ES (Z5 fo-of ss, uu ,z,t)dzdodt — c f; e 

Of E T(R, pz t)dzdpdt + ——- —_ x f° - mye T,(0, b, z, t)dzdpdt -—_ f° (0 

t) p f tQ. iz, t)[æzc + Cs CR, , Z,t)] — araic} X [v (5, 6,7, t) - 

9, (R, d, z, O]dzdédt — fj (0 —0 f; Jy" Fille + C (rd zt] 7 arac): [ns zt) — 
B,(r, ,z,t)] X T (r, dz, tO dzdérdrdt —— — =z f"(8 — t) fe r f7 [ri o, L, t)(aac + 


Co) — aras c] dédrdt (fo (0 — t) x fo" f^ Dv (0, d, z t) — 9,0, 4, z, ] Zilatac + 
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O pR 2 L = 
GR, bz, €] dzdodt — +1- fo for fo" f. Qd zi 0 — 9. o, z t) Eile + 
Cis(R, ,2, t] dzdddr(0 — t)dt a + 20( azic + Co) 2). eic = Saaz yO - 


D fr J, [Card zt) — Cur p, 2,0) |dgdrdt. 





3. RESULTS AND DISCUSSION 

Now we analyzed redistribution of reagents and temperature in the considered reactor during growth 
of films to formulate conditions for improvement of properties of films. Figure 2 illustrates dependence of 
concentrations of reagents on substrate's rotation frequency. Curve 1 shows such dependence at atmospheric 
pressure without accounting natural convection. Curve 2 shows such dependence of concentrations of 
reagents on substrate’s rotation frequency at decreased in 10 times pressure. Curve 3 describes the 
dependence of the concentration of reagents on substrate’s rotation frequency at atmospheric pressure with 
account natural convection. The figure shows, that increasing of rotation frequency leads to increasing of 
homogeneity of epitaxial layer. Figure 3 illustrates dependence of concentration of reagents on diffusion 
coefficient D in gas-carrier. Curve 1 shows dependence of the concentration of reagents at atmospheric 
pressure without accounting natural convection. Curve 2 shows dependence of the concentration of reagents 
at decreased in 10 times pressure without accounting natural convection. Curve 3 shows dependence of the 
concentration of reagents at atmospheric pressure with accounting natural convection. This figure describes 
monotonous decreasing of considered concentration. Figure 4 illustrates dependence of concentration of 
kinematic viscosity of gas-carrier without accounting natural convection. Curves, analogous corves 2 and 3 
on previous figures, have small differences with the presented curve. Figure 5 illustrates dependence of 
concentration of reagents on velocity of mixture of reagents and gas-carrier at the entrance to the reaction 
zone Vo. Curve 1 shows dependence of the concentration of reagents at atmospheric pressure without 
accounting natural convection. Curve 2 shows dependence of the concentration of reagents at decreased in 10 
times pressure without accounting natural convection. Curve 3 shows dependence of the concentration of 
reagents at atmospheric pressure with accounting natural convection. This figure describes monotonous 
decreasing of considered concentration. 
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Figure 2. The dependence of the concentration of the Figure 3. The dependence of the concentration of the 
mixture of gases on the frequency of rotation of mixture of gases from its diffusion coefficient 
the substrate 
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Figure 4. The dependence of the concentration of the mixture of gases from the kinematic viscosity 
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Figure 5. The dependence of the concentration of the mixture of gases from its input velocity 


It is attracted an interest choosing of power of induction heating of reaction zone. This choosing 
should compensate losses of heat due to convective heat transfer [21]. In this situation relaxation time of 
temperature could be estimate by using recently introduced criterion [22]. This time is equal to 9 «(67 -1) 
R?/24%, where 2 is the average value of the thermal conductivity coefficient. In this case, the power 
required to compensate for the cooling of the region of the formation of the epitaxial layer can be estimated 
from the following relation 


fo r:pG,d.z Odr = o: T*(R,ó,z, t) + O-v,(R, Quz, O /AnLR? [16]. 


4. CONCLUSION 

In this paper based on considered analytical approach for analysis of mass and heat transfer during 
film growth in reactors for gas phase epitaxy we analyzed this growth with account natural convection and 
chemical interaction between reagents. In this situation we analytically analyzed growth of films from gas 
phase in more common case in comparison with cited references to increase predictability of this 
technological process (it should be noted, that analytical approach for modeling are usually more 
demonstrative in comparison with numerical one and leads to decrease calculation time). As a result of this 
analysis, the obtain conditions to increase homogeneity of the grown epitaxial layers with changing of values 
of parameters of growth process with higher exactness. 


Int. J. of Adv. in Appl. Sci. Vol. 8, No. 3, September 2019: 195 — 203 


Int. J. of Adv. in Appl. Sci. ISSN: 2252-8814 o 203 





REFERENCES 

1] I. P. Stepanenko, Basis of microelectronics. Moscow: Soviet radio, 1980. 

2] V.G. Gusev and Yu M. Gusev, Electronics. Moscow: Higher School, 1991. 

3| V.I. Lachin and N. S. Savelov, Electronics. Rostov-on-Don: Phenics, 2001. 

4] A. A. Vorob'ev, et al., “The use of magnesium to dope gallium nitride obtained by molecular-beam epitaxy from 
activated nitrogen," Semiconductors, vol. 37(7), pp. 838-842, 2003. 

5] L. M. Sorokin, et al., “Electron-microscopic investigation of a SiC/Si (111) structure obtained by solid phase 
epitaxy," Technical Physics Letters, vol. 34(11), pp. 992-994, 2008. 

6] V. V. Lundin, et al., “Effect of carrier gas and doping profile on the surface morphology of MOVPE grown heavily 
doped GaN: Mg layers," Semiconductors, vol. 43(7), pp. 963-967, 2009. 

7] Y. E. Bravo-García, et al., “Growth and characterization of InAsSb layers on GaSb substrates by liquid phase 
epitaxy," Mat. Sci. Sem. Proc., vol. 40, pp. 253-256, 2015. 

8] Y. Li, L. E. Antonuk, et al., “Effects of x-ray irradiation on polycrystalline silicon, thin-film transistors," J. Appl. 
Phys., vol. 99(6), pp. 064501, 2006. 

9] A. Chakraborty, et al., *Nonpolar [l-plane p-type GaN and p-n-junction diodes," J. Appl. Phys., vol. 96(8), 
pp. 4494, 2004. 

10] H. Taguchia, et al., “Evaluation of crystallinity of GaN epitaxial layer after wafer dicing,” Mat. Sci. Sem. Proc., 
vol. 41, pp. 89-91, 2016. 

11] M. Mitsuhara, et al., “Beryllium doping of InP during metalorganic molecular beam epitaxy using 
bismethylcyclopentadienyl-beryllium," J. Cryst. Growth, vol. 183(1-2), pp. 38-42, 1998. 

12] R. A. Talalaev, et al., “On low temperature kinetic effects in metal-organic vapor phase epitaxy of III-V 
compounds," J. Cryst. Growth, vol. 230(1-2), pp. 232-238, 2001. 

13] H. S. Carslaw and J. C. Jaeger, Conduction of heat in solids. Oxford: At the Clarendon Press, 1964. 

14] Yu. D. Sokolov, *About the definition of dynamic forces in the mine lifting," Applied Mechanics, vol. 1(1), 
pp. 23-35, 1955. 

15] E. L. Pankratov and E. A. Bulaeva, “On optimization of technological process to decrease dimensions of transistors 
with several sources," Micro and Nanosystems, vol. 8(1), pp. 52-64, 2016. 

16] E. L. Pankratov and E. A. Bulaeva, “An analytical approach for analysis and optimization of formation of field- 
effect heterotransistors," Multidiscipline Modeling in Materials and Structures, vol. 12(4), pp. 578-604, 2016. 

17] E. L. Pankratov, “On optimization of manufacturing of power amplifier circuit based on bipolar heterostructures to 
increase density of their elements. Influence of miss-match induced stress," Advanced science, engineering and 
medicine, vol. 9(10), pp. 849-863, 2017. 

18] E. L. Pankratov, “On decreasing of dimensions of transistors with two sources by optimization of technological 
process," Journal of computational and theoretical nanoscience, vol. 14(10), pp. 4947-4954, 2017. 

19] E. L. Pankratov and E. A. Bulaeva, “Doping of materials during manufacture p-n-junctions and bipolar transistors. 
Analytical approaches to model technological approaches and ways of optimization of distributions of dopants," 
Reviews in Theoretical Science, vol. 1(1), pp. 58-82, 2013. 

20] G. Korn and T. Korn, Mathematical Handbook for scientists and engineers. Definitions, theorems and formulas for 
reference and review. New York, Second edition. McGraw-Hill Book Company, 1968. 

21] E. L. Pankratov and E. A. Bulaeva, “On prognosis of epitaxy from gas phase process for improvement of properties 
of films," 3D research, vol. 6(4), paper 40, 2015. 

22] E. L. Pankratov and E. A. Bulaeva, “Optimal criteria to estimate temporal characteristics of diffusion process in a 
media with inhomogenous and nonstationary parameters. Analysis of influence of variation of diffusion coefficient 
on values of time characteristics," Reviews in Theoretical Science, vol. 1(3), pp. 305-316, 2013. 








On approach for homogeneity increasing of films grown from the gas phase with ... (E. L. Pankratov) 


